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Abstract 

An exact statistical mechanical derivation is given of the critical Casimir forces for Ising strips 
with arbitrary surface fields applied to edges. Our results show that the strength as well as the 
sign of the force can be controled by varying the temperature or the fields. An interpretation of 
the results is given in terms of a linked cluster expansion. This suggests a systematic approach for 
deriving the critical Casimir force which can be used in more general models. 
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Casimir forces [l| arise in the quantum electrodynamics of geometrically restricted sys- 
tems, for instance between two metal plates in vacuo, because the photon spectrum is 
modified; typically the force is attractive. Fisher and de Gennes 3] proposed that anal- 
ogous Casimir forces should arise in condensed matter systems near a second-order phase 
transition, the agent being thermally-excited fluctuations of the order parameter, e.g. of 
the density, rather then of the photon field. Of particular interest, both experimental 
and theoretical, is their scaling-theoretic prediction that such interactions should have a 
power law dependence on distance in the critical scaling region. For example, in spa- 
tial dimension d for walls separated by a distance N, the Casimir force per unit area is 
Fcas — ~N~ d, &(N/^), where £ is the bulk correlation length j^, ^, 4]. (All free energies and 
forces are expressed in units of /cgT.) The tunability of critical Casimir interactions will be 
crucial for many applications in micro and nano systems, e.g. in colloids or in various micro 
or nano-electromechanical devices, in particular, to be able to produce repulsive interactions 
to counteract the omnipresent attractive Casimir quantum electrodynamical force. 

If the system is confined to a film, one would expect on intuitive grounds that the ge- 
ometrical effect on the order parameter fluctuations of the strip boundaries, would be to 
reduce the entropy; thereby establishing a strictly repulsive force. This argument is incom- 
plete because some energetic factors are neglected. We propose an alternative view that 
admits attractive forces. Let us examine the following approximate treatment of a 2d Ising 
strip of finite width N and with free boundaries. Following Privman and Fisher 5|, the 
low energy excitation of such a model are domain walls, each with free energy Nt where 
t is the incremental free energy per unit length for an interface perpendicular to the strip 
axis. A collection of these is then treated as a Id Ising model with cyclic boundary condi- 
tions. A simple calculation shows that the limiting incremental free energy, per unit length 
of strip, is f x = — log (l + e~ Nr ) and then the Casimir force, again per unit length, is 
T Cas = -df x /dN = -r/(e Nr + 1). In the scaling limit Nt = x, N -> oo, r -> 0, (CQ) 
becomes 

NFcasi*) - (1) 

LI 

Notice that this force is attractive, but the power of iV disagrees with [2J . What is missing 
in this formula? Firstly, there is no mention of any thermally-excited intrinsic structure 
of the individual domain walls. This could perhaps be included for the planar Ising model 
with free boundaries by calculating the partition function of a single interface connecting 
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FIG. 1: (color online) The scaling function y±, y 2 ) of the critical Casimir force (J7]) for the 

isotropic lattice with K\ = K 2 , y\ = h\N = 00 and for the several values of y<i = h 2 N. 

the edges with one end fixed, but the other free; this merely reproduces §1$). Secondly, the 
detailed interactions of these domain walls must be quite complicated, but they are evidently 
accounted for in exact calculations on planar Ising strips. The first known results ||]] were 
for the strip with zero bulk field and either free boundary condition or fixed boundary spins, 
both the ++ and H — conditions. 

The partition function for the free strip is [3] 

Z f = [cosh Nj{uj) + sinh N-f(co) cos S'(u)} . (2) 

0<u;<7r 
e iNu =-l 

The Onsager functions [s| 7(0;) and 5' (u) are given by cosh7(c<j) = cosh 2K 2 cosh 2K* — 
sinh 2K 2 sinh 2K\ cos u, and 

where z = e luJ , A = exp 2(Ki + K 2 ) and B = exp 2(Ki — K 2 ). K x and K 2 are the nearest 
neighbour couplings (in units of ksT) in (0, 1) and (1,0) direction, respectively. K* is the 
dual coupling given by the involution sinh(2J^) sinh(2J^*) = 1. The bulk free energy is given 
by extracting a factor e Nl /2 from the argument of the logarithm leaving an incremental free 
energy. Removing an additional ^-independent factor gives 

/ x = f^\n[l + (l/d'(u)e-™^], (4) 
where d'(u) = (1 + cos 5' (lj)) / (1 — cos5'(u)). The Casimir force per unit length is now 

?Ca.(N,T) = -/^ KMe 7 2 ^l ) + 1] - (5) 



Comparing this with ([T]), we see that ([5]) also describes an attractive force, but there is also 
an integration, r is replaced by 7(0;) under the integral sign and we have exp 2N r y(u) rather 
than expiVr in the denominator. In addition, there is the prefactor d'(u) multiplying this 
term. Taking the scaling limit of (jSJ) then recaptures the Fisher-de Gennes law with the 
correct power. 

The aims of this Letter are two-fold: firstly, we show how (jSj) is generalised to include 
non-zero surface fields, which allows us to control the sign of the Casimir force at will. 
Secondly, we interpret this result as a linked cluster expansion and then indicate why it 
might be appropriate for models more general than the planar Ising ferromagnet with nearest 
neighbour interactions. 

A fortnalisatiou of the technic of Schult, Mattis and Lieb aUows us to calculate 
the partition function of a cylindrical lattice with circumference M, height N, with its axis 
in (0, 1) direction, and end fields hj > 0, j = 1,2 in a straightforward way. The fields are 
introduced by taking a free-edged cyclic strip and adding an extra ring of spins at each end; 
these spins are forced to take the value +1. These fixed spins are then coupled to the free 
lattice by bonds of strength hi at the bottom and h 2 at the top. A different technique will be 
needed if h±h 2 < 0, as will be seen. The Casimir force per unit length in the (1, 0) direction 
as M —* 00 is obtained from the incremental free energy by taking the derivative in respect 
to N: 

Fcasihi, h 2 , N,T) = - [ W — (6) 

where d*(u) = (1 + cos<J*(u;))/(l - cos<5*(u;)) and A s = (e" 7 ^ - wj) / (e^ - w 3 ) , j = 
1,2. e l5 *^ is given by (J3j) with B replaced by B~ x . The values W\ and w 2 which are the 



wetting parameters for the force are given by 



10| Wj = e 2K2 (cosh2K! - cosh 2^) /sinli2iif 1 . 



It is crucial to note that Aj(u) can take both positive and negative values; this is why either 
sign of the Casimir force is possible in principle. In the scaling limit, N — > 00,7(0) = K 2 — 
K{ -> such that x = Nj(0)sgn(T-T c ) is fixed (as t = (T-T c )/T c -> 0, K 2 -K\ ~ -AK c t) 
and Nu = u leads to Tcas{h\, h 2 , N, T) = N~ 2 $ + (r), with 

Jn X LQj A i X 2 r 2\(x.u) 1 1 
X+(0)X^X2 

where r = (x, y±, y 2 ) with yj = h 2 N, X(x, u) = \/x 2 + u 2 and X ± (y) = X(x, u)^f(x — 2e 2Kc y), 
Xf = X ± (yj),j = 1, 2. For T < T c , 7(0) is the surface tension in the (0, 1) direction of the 
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FIG. 2: (color online) The scaling function y±, y 2 ) of the critical Casimir force (J7]) for the 

isotropic lattice with K\ = K 2 y\ = h\ N = and for the several values of y 2 = h\N. 

2d Ising model; it is the inverse correlation length for T > T c . At x = ([7]) reduces to the 
universal Casimir amplitude, which equals — 7r/48 for both hj = and hj = 00. Interesting 
examples of the scaling function $+(r), which demonstrate that the critical Casimir forces 
can switch from attration to repulsion by varying the temperature, are shown in Figs. [1] and 
[2j They were evaluated numerically from (J7|) for several choices of the scaling variables y± t2 ■ 
The case with h\h 2 < can be approached from that with h\h 2 > by reversing the end 
spins between x — 1 and x = s + 1 on one face of the cylinder thereby creating an interface 
with terminations in the same face. This is followed by taking the limit as M — > 00, as 
before. With hj > 0,j — 1,2, we find the ratio of partition functions for strips with and 
without the interface to be: 

_ s/x r duj ie lsw ttm(5*/2) (g+ - B~ A 2 e~ 2N ^) 
6 ~ J T 2ti l + A 1 A 2 e~ 2N ^ ( ) 

where B ± = [e^ 1 ^ — e~ iK2 Wij /(e 7 ^-* — wi). We are interested in the limit s — >• 00 of 
the rhs of (jSJ) per unit length. The asymptotics for large s is dominated by the nearest 
singularity to the real axis in the strip — tx < uj < it. The branch cuts associated with 
sinh7(u;) do not occur and poles are simple zeros of the denominator of (jSJ). Fortunately, 
the problem can be related to the diagonalisation problem of the transfer matrix in the 
direction (1, 0) UJ (here we are transferring in (0, 1) direction) by looking for the solution 
in the variable k such that uj = ij(k) and lim e ^ + ± e) = ±ik, where the function 

j(k) is defined as the Onsager function, but with K\ and K 2 interchanged. Then finding 
zeros of the denominator of (jSJ) becames equivalent to solving the spectrum discretisation 
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condition for the strip transfer matrix in the (1,0) direction, which was studied in detail in 



Ref. 



Ill: 



e 2i(N+l)k = c 2^(fc) e *y- leifew 2-l , Q s 
e ik _ Wl e ik _ W2 

e %8 (k) j g obtained from e lS '^ by interchanging K 1 and K 2 . In the scaling limit N — > oo, 
K 2 — K* — > and N(K 2 — K*) = x, we find f x = (l/N)X(x, u ). Hence the solution for the 
Casimir scaling function i?_(r) = -# x (r) + $+(r) has the implicit form with 

< io) 

N 2 X(x, u ) 

where Mo(zi) solves the quantisation condition in the scaling limit 

e 2 » = fin 

z-(o)z+z+' 1 J 

where Zj* 1 = Z ± (yj),j = 1,2 is derived from -X±(?/) by replacing A(x,m) by zu. The deriva- 
tives of no can be calculated straightforwardly from ffTTj) . In Figs. [3] and H] we plot t? x as a 
function of x evaluated numerically for some choices of the scaling variables y\ and y 2 . Our 
results for the special case of yi = y 2 agree with those reported in Ref. the change of 
sign of the scaling function i9 x is associated with the localisation-delocalisation transition 



13j| . This feature remains for a slightly broken symmetry, i.e., for y\ y 2 and y\^ 2 small. 
For strongly asymmetric strips the excess scaling function of the critical Casimir force is 
always positive. 

We now interpret (jSJ) in terms of statistical mechanical ideas. Expanding the integrand 
gives 

where Cj(u) = ^ tan 5 ^ , j = 1,2. Although this is not immediately apparent, this is in 



fact a linked cluster expansion as we now show. Equation ([12]) can be understood by going 
back to the partition function formula in terms of a transfer matrix V 

Z= (b 1 \V N \ b 2 ) (13) 

where | bj) describes the edge state with field hj,j = 1,2. Instead of using the Schultz, 
Mattis and Lieb technology, this can be developed by expanding with basis of eigenvectors 
of V giving 

Z ~ e -*E£i7("j) 

IF = EE f2 „ V ^ I ("hn){(")2n I h), (14) 
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where | (u>) 2n ) denotes a 2n-fermion eigenstate of V 9fl. This would certainly not be the 
chosen way of obtaining ([6]), since we would need to evaluate the matrix elements (bi | 
0^)2n)> (( UJ )2n I b 2 ); this has been done with some effort and the result is typically Wick- 
theoretic in form: 

2n 

(b, I (u) 2n ) = Y,(- l ) m fM^m)(b 3 I A lm (o;) 2n ) (15) 

m=2 

where Ai m (u;)2 n = (uj 2 , ■ ■ ■ w m -i, <*Wi ■ ■ ■ u 2n ) and fj{u>i, uj 2 ) = iCj^S^-^ is tne contrac- 
tion function or, alternatively, a scattering matrix element for a pair of fermions off the 
wall described by | bj). Notice that since the Aj (oj) are even, but tan|<5*(u;) is odd, the 
contraction is antisymmetric as it should be for fermions. Thus, ffl5|) is a Pfaffian [l4]. The 



graphical representation of (TH1) and (Tl5l) is discussed in 15|. For each n we have a weighted 
sum of disjoint loops, each having an even number of vertices, the vertex weight Cj(u) and 
the Kronecker delta edge weight. The occurrence of the Kronecker delta in the contraction 
function is mandated by translational symmetry. Thus, the multiple sum for each loop 
becomes just a single sum on implementing the deltas. Asymptotically as M — > oo, each 
such sum is to leading order M times a single integral. We can now apply the linked cluster 
theorem to exponentiate ( fl4|) . Eq. ( fl2l) is recaptured, for the excess free energy per unit 
length in the (1, 0) direction since the factor of 1/n in (1121) comes directly from a symmetry 
number argument [if]]. Each term is then to be thought of as a weight of a "loop" with 
2n vertices. The loop is reflected n times off the upper boundary and n times off the lower 
boundary, with "momentum" conservation at each reflection; thus n maybe thought of as a 
topological quantum number. Starting from ( fl3l) . ( fl4l) and ( fl5|) . we have re-derived (fl2l) . in 
a way which allows us to identify the multiplier of exp 2^7(0;) in (jHJ) as a product of two 
scattering matrix elements, one from each edge. Clearly 7(0;) in ( Tl4|) is a Fermion energy. 
Thus we have a complete intuitive understanding of ([6]). We can take the scaling limit ei- 
ther in (fl2|) or ([6]) (as we have already done) with the same outcome. This procedure even 
converges after taking T — > T c in either (fT2"l) or ([6]), since then 7(0;) oc| u \. 

Two approximation schemes are in order. Firstly, we could consider how well partial sums 
of the virial series ( Tl2l) approximate the exact result, so that we can assess the contribution 
of the different reflection number sectors to the result. Secondly, in the scaling limit the one- 
particle energy should be universal. The same is not likely to be true for the scattering matrix 



elements. Correlation droplet theory 



171 ] provides an approximate method for calculating 



FIG. 3: (color online) The excess scaling function i9 x (x, 2/1,2/2) of the critical Casimir force f) 1 Of) 
for the isotropic lattice with K\ = K2 y\ = h\N = 00 and for the several values of 7/2 = h\N . 

them and thus for extending the scope of our results. 

In this Letter, we have described new, general, exact results for the critical Casimir force 
in a planar, rectangular Ising ferromagnet with applied fields hi and h 2 on the edges. Each 
field can have arbitrary sign and magnitude. Both with h\h2 > and with hiti2 < 0, we 
show that the force can be attractive or repulsive, according to the tuning of the parameters. 
The compensation of attractive, quantum van der Waals forces which this will allow has 
implications which may well prove crucial for applications. Mean field calculations are in 



qualitative agreement with our results 

n 



181 ] . There are also related results from the continuum 
model We also interpret the representation of the Casimir force as in ([5]), ([6]) and (fl2|) 
in terms of the linked cluster expansion. This suggest an associated droplet picture which 
enhances the original finite size scaling ideas of Privman and Fisher j^j in this context; this 
will also give new, systematic approximations for calculating critical Casimir forces in planar 
systems and perhaps even in d — 3. Our results can be directly applied to e.g. 2c? binary 
fluid membranes with protein inclusions close to the demixing point {2^. 
DBA acknowledges Max-Planck-Gesellschaft for hospitality. 
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